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The thermal entanglement of a two-qubit anisotropic Heisenberg XY Z chain under an inhomoge- 
neous magnetic field b is studied. It is shown that when inhomogeneity is increased to certain value, 
the entanglement can exhibit a larger revival than that of less values of b. The property is both 
true for zero temperature and a finite temperature. The results also show that the entanglement 
and critical temperature can be increased by increasing inhomogeneous exteral magnetic field. 
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I. INTRODUCTION 



II. THEORETICAL TREATMENT AND 
RESULTS 



As a valuable resource in quantum information and 
quantum computation [1] , quantum entanglement has at- 
tracted numerous attention over past decade years. In or- 
der to realize the quantum information process, a great 
deal of effort has been devoted to study and characterize 
the entanglement in solid state systems [2,3]. A typical 
example in solid system is the spin chains which are the 
natural candidates for the realization of the entanglement 
compared with the other physics systems. The Heiser- 
berg model can describe interaction of qubits not only 
in solid physical systems but also in many other systems 
such as quantum dots [4] and nuclear spin [5]; therefore, 
there are numerous studies on the Heisenberg models [6- 
13]. It is turned out that the critical magnetic field B c is 
decreased with the increase of the anisotropic parameter 
7 but the critical temperature T c is improved. When B 
crosses B c the concurrence C drops suddenly and then 
undergoes a "revival" for sufficiently large 7 in Ref [14], 
but it only discussed the uniform magnetic field case. As 
we all know, in any solid state construction of qubits, 
there is always the possibility of inhomogeneous Zeeman 
coupling [15,16]. Quite recently the effect of inhomoge- 
neous magnetic field on the thermal entanglement of an 
isotropic two-qubit XXX spin system has been studied 
in Ref [17]. Including interaction of Z-component J z , 
Ref [18] also studied the effect of inhomogeneous mag- 
netic field on the thermal entanglement in a two-qubit 
Heisenberg XX Z spin chain. 



However, the entanglement for a XYZ spin model un- 
der an inhomogeneous magnetic field has not been dis- 
cussed. Therefore, in this paper we investigate the influ- 
ence of an external inhomogeneous magnetic field on the 
entanglement of a two-qubit Heisenberg XYZ system at 
thermal equilibrium. Our studies show that the inhomo- 
generous extral magnetic field can make the revival en- 
tanglement larger, improve the critical temperature, and 
enhance entanglement. 



The Heisenberg Hamiltonian of a N-qubit anisotropic 
Heisenberg XYZ model under an inhomogeneous mag- 
netic field is 



H 



1 N 

i=l 



+(B + b)a* i +(B-b)al +1 ], 



(1) 



where (erf , g\, erf) are the vector of Pauli matrices and 
Jj (j = x, y, z) is the anisotropic coupling coefficient be- 
tween the nearest two spins. The parameter Jj > 
means that the chain is antifcrromagnetic, and ferromag- 
netic for Jj < 0. The magnetic fields on the nearest- 
neighbor two qubits are B — b and B + b, respectively, 
the value of b controls the degree of inhomogeneity. 

For a spin system in equilibrium at temperature T, 
the density matrix is p — (1/Z) exp(— H/keT), where 
H is the Hamiltonian of this system, Z is the partition 
function and ks is the Boltzmann constant. Usually we 
write fee = 1' For a two-qubit system the thermal en- 
tanglement can be measured by the concurrence C which 
can be calculated with the help of Wootters' formula[l9j 

C = max(0, 2maxAi — Y^t=i ^»)' wnere A» is the square 
roots of the eigenvalues of the matrix 



R = p{a\®al)p*{al®al), 



(2) 



where the asterisk indicates complex conjugation, the 
concurrence C ranges from zero to one. 

Consider now the Hamiltonian H for the anisotropic 
two-qubit Heisenberg XY Z chain in an inhomogeneous 
magnetic field. The Hamiltonian can be shown as 

H = J(cr+(7^ + a^cr+) + Jj(a+a+ + +a^a^) + y' 7 !' 7 
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where J = ^+ Jy) , 7 



and a ± = \{a x ± ia v ). 
Among these parameters er* are raising and lowering 



J x J y 
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operators respectively and 7 (0 < 7 < 1) measures 
the anisotropy in the XY plane. In the standard basis 
{|00), 1 01) , 1 10), |11)}, the Hamiltonian can be expressed 



as 
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The eigenvectors and eigenvalues of H are easily ob- 
tained as following forms H\ip ± ) = (— 4f ± OlVw an d 
i/|E ± ) = (4f ± ?7)|S ± ), with the eigenstates are |V' ± ) = 
jV±[Mfl|01) + |10>] and |£±) = M±[^±^|00) + |11)], 
where r] = \J B 2 + J 2 j 2 and £ = V& 2 + J 2 , the nor- 
malization constants are TV* = \j \J\ + (b ± £) 2 / J 2 and 
M ± = l/y/l + (B ± 7y) 2 /J 2 7 2 . One can notice that 
these four eigenstates are all entanglement states when 
J 7^ so it means that entanglement exists for both an- 
tiferromagnetic (J > 0) and ferromagnetic (J < 0) cases. 

Now we can calculate the thermal entanglement of this 
system. The density matrix can be written as the follow- 
ing form 



P = 



/M+ v \ 

wi z 

z uj 2 

V u n_ J 



The square roots of the eigenvalues of the matrix R are 
Ai,2 = IVM+M- ± v\, A 3i4 = |V w i^2 ± z|- (4) 

The exact values of these nonzero matrix elements can 
be obtained by knowing the spectrum of H. We obtain 

/j + = ^-e~~2~P[cosh(r)/3) — — sinh(r//?)], 



/i_ = * /3 [cosh(r ? /3) + -sinh(r//3)], 
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* = -^-e^sinh(£/3), 



Wl = i e ^[cosh(^)-|sinh(^)], 



0J2 = |e* /3 [cosh(£/3) + I sinhfoS)], 



u = -— e '2 s ^ sinh(r//3), 



(5) 



where Z is the partition function given by 

Z = tre- fm = 2[e-^ 13 cosh(7?/3) + e^ 13 cosh(£/3)]. (6) 
Thus from Eq.(4)-(6) we find that 



A1.2 = ^e-^Jl + ^ sinh^/?) 2 T ^ sinh^)!, 



\3,4 




^sinh(C/3) 2 T|sinh(e/3)|. (7) 



We can calculate the thermal entanglement which is mea- 
sured by the definition of the concurrence. One can 
notice that Eq.(7) arc the same as the case of b = 
[20]. The concurrence which is directly determined by 
Xi(i = 1,2,3,4) is invariant under the substitutions 
7 — ► —7 and J — ► — J but the eigenvalue is variant 
with the substitution J z — ► — J z so that the concur- 
rence is also different. We will take the cases J > and 
< 7 < 1 into account in this paper. 

From Eq.(7) and the definition of the concurrence, we 
derive the C for T = as following 



C(T = 0) = 



£,<v-Jz 



(8) 



Although the parameters J, 7, 77 and £ are independent 
of J z in the case of two interacting qubits, the value of 
J z is very important in determining the point of the 
piecewise function so that it can play role in the pairwise 
entanglement. We will show it in Fig.l. 

The concurrence C as a function of b at T — for 
three values of 7 are given in Fig.l. Consider first the 
nonuniform magnetic field (B = 0) case [Fig. 1(a)]. With 
increasing b, the concurrence C is initially constant and 
equal to its maximal value 1. It then drops suddenly 
as a critical value b c is reached, where the critical in- 
homogene ous magnetic fi eld b c is given by £ = rj — J z 
or b c = ,J(ri - J z ) 2 - J 2 . At the critical point (T = 0, 
b = b c ), the entanglement becomes a nonanalytic func- 
tion of b and a quantum phase transition occurs, but for 
b > b c the concurrence C undergoes a revival before de- 
creasing to zero. The critical inhomogeneous magnetic 
field b c is increased by increasing the anisotropy param- 
eter 7, which means the region that entanglement keeps 
its maximum value is broaden with the increasing 7. 
However, for fi / case, entanglement exhibits novel 
property shown in Fig. 1(b). For certain relationship of 
the parameters, the value of revival can be larger than 
its value before dropping, for example, Figl.(b) 7 = 0.2 
curve. In addition, the region of entanglement keeping 
constant and entanglement values is increased with the 
increasing 7. In Fig.l (c), we further show the role of J z 
in existing larger revival. One can see clearly the values 
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FIG. 1: The concurrence in the X Y Z model is plotted vs b for 
various values of anisotropy parameter. (a):From left to right 
7 = 0.2 , 0.6 ,0.9, respectively, where B = and J z = — 1. 
(b): From bottom to up 7 = 0.2, 0.6, 0.9, respectively, B — 0.8 
and J z = — 1. (c) gives the concurrence with different values 
of J z : J z — —0.2 (solid line), J z = —0.6 (dashed line) with 
7 = 0.6 and B = 0.8. For all plotted T = and J = 1. 




FIG. 2: (color online) Concurrence versus the inhomogeneous 
magnetic field field 6 and J z , where J = 1, B = 4, T = 0.2 
and 7 = 0.3. 



of J z broaden the region of entanglement keeping con- 
stant values. The most important is that only when J z 
achieve certain values, the value of entanglement revival 
can larger than that of before dropping. 

Now we discuss the conditions of existing revival and 
larger revival. From Eq.(8), we see if £ ^ 77 — J z , i.e., 
b c ^ vCT^ Jz) 2 — J 2 , the concurrence C will exhibit 
revival phenomenon. But if we need larger revival value, 
combining j > ^ and £ ^ T] — J z , we have J z > 77 — 2 . 
Therefore, the condition existing larger revival is as 



> V(v - J z ) 2 - J 2 



(9) 



In figure 2 we give the plot of concurrence as a func- 
tion of b and J z at a fixed temperature and magnetic field 
for 7 = 0.3. It is shown that the critical inhomogeneous 
magnetic field b c is increased with the decrease of the in- 
teraction of z-component J z and there is no entanglement 
at this critical point. As b increases, the concurrence C 
is initially constant which is nearly invariant by chang- 
ing the interaction of z-component J z . It then decreases 
suddenly as the critical inhomogeneous magnetic field b c 
is reached. However, instead of vanishing for b > b c , 
the concurrence persists and undergoes a revival before 
decreasing. The most important thing is that the value 
of entanglement revival can larger than that of before 
dropping. The maximal entanglement value exists in the 
revival region, which means we still can observe larger 
revival phenomenon for thermal entanglement. 

With B = 4, the concurrence as a function of b and T 
are given in Fig. 3. If b is lower than a certain value, there 
are two areas showing entanglement, and the entangle- 
ment first decreases and then undergoes a revival before 
decreasing to zero. But when b is larger than this certain 
value, there is no revival phenomenon and the entangle- 
ment is decreased monotonously with the increase of T. 
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FIG. 3: (color online) Concurrence versus temperature and 
inhomogeneity for 7 = 0.2 and B — 4. The coupling constant 
J — 1 and J z = 1. 
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FIG. 4: (color online) The concurrence in the XYZ spin 
model is plotted vs b and anisotropy parameter 7, where 
T = 0.4 and B = 0.8. The coupling constant J z = -0.6 
and J — 1. 



However, in the region of larger (b, T) parameter space, 
we can derive larger entanglement. In addition, a crit- 
ical temperature T c is improved with the increase of b. 
The property opens out the role of inhomogeneous mag- 
netic field b in improving entanglement. Ref.[10] pointed 
out that a pairwise entanglement in iV-qubit isotropic 
Heisenberg system in certain degree can be increased by 
introducing external field B. Our study shows the inho- 
mogeneous magnetic field b can have a better action on 
improving entanglement for thermal entanglement. 

Figure 4 shows the entanglement as measured by 
the concurrence in terms of the inhomogeneity b and 
anisotropy parameter 7 at a finite temperature T for a 
finite external magnetic field B (B = 0.8). We can ob- 
serve that entanglement exists in two regions: the one is 
in large 7 and small 6; the other is in small 7 and b. If 
the anisotropy of XY plane 7 is large enough, we just 
need small inhomogeneous external field b or even with- 
out inhomogeneity, vice versa. At this point, the role of 
b is something similar to the anisotropy XY plane. It 
seems to be possible to control the anisotropy by adjust- 
ing external magnetic field. 



III. CONCLUSION 

In conclusion, we have studied the thermal entangle- 
ment in an anisotropic two qubits Heisenberg XYZ chain 
under an inhomogeneous magnetic field. Through ana- 
lyzing the T = case, we find that conditions of existing 
revival phenomenon and larger revival. The larger re- 
vival phenomenon still exists at fixed finite temperature. 
If the parameters are proper and the inhomogeneity is 
large enough, in the revival region the entanglement is 
larger than that of without inhomogeneous magnetic field 
for fixed temperature. Then we show by increasing the 
inhomogeneous magnetic field b the entanglement and 
the critical temperature T c can be improved. 
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